Opacity and Entanglement of Polymer Chains 
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We argue that the mean crossing number of a random polymer configuration is simply a measure 
of opacity, without being closely related to entanglement as claimed by several authors. We present 
an easy way of estimating its asymptotic behaviour numerically. These estimates agree for random 
walks (theta polymers) , self-avoiding walks, and for compact globules with analytic estimates giving 
log N,a — b/N 2 " -1 , and N 1 ' 3 , respectively, for the average number of crossings per monomer in the 
limit N — > oo. While the result for compact globules agrees with a rigorous previous estimate, the 
result for SAWs disagrees with previous numerical estimates. 
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Topological properties of linear polymer chains are, 
strictly spoken, only defined for closed rings where the 
theory of knots can be applied with several non-trivial 
results But even for open chains quantities like 

linking, twist and writhe can be applied. Important ap- 
plications are supercoiling of DNA, DNA electrophoresis 
H , and the supposed absence of knots in protein back- 
bones Q . Another class of problems where entanglement 
should play a crucial (but so far not yet fully understood) 
role is the rheology of semidilute solutions. 

One important quantity for studying entanglement is 
the writhe which is defined as the number if signed cross- 
ings in a projection of a 3-d non-selfintersecting oriented 
curve. If two parts of the curve seem to cross when seen 
from a particular angle, this crossing contributes +1 or 
— 1 to the writhe, depending on whether the direction of 
the front part is obtained by a right or left turn from the 
direction of the part behind it. Its interest stems from 
the fact that for closed loops it is related to linking which 
is a topological invariant [||. 

The number of crossings C was introduced in (^^J as 
a simplified version of the writhe. In it, all crossings con- 
tribute with the same sign. Of particular interest is its 
average value (C), averaged over all angles of projection, 
called the mean crossing number. In these papers it was 
shown that for a self avoiding random walk (SAW) of N 
straight bonds 



(C) - N a for N 



(1) 



with 1 < a < 2. Numerical simulations gave a = 
1.122 ± 0.005, but it was argued that this might actu- 
ally be a lower estimate, the true value being higher [Q. 

In later simulations, Arteca [|| found a value 1.20 ± 
0.04 for SAWs and 1.34 to 1.4 for protein backbones [Q. 
Indeed an increase of the value of (C) had also been seen 
in |lj for SAWs with self-attraction, and it was conjec- 
tured in [[jJD that (C) is a useful observable for detecting 
the coil-globule transition. Due to its supposed impor- 
tance, (C) was called the "entanglement complexity" in 
H , and was shown there (by non-rigorous arguments) to 



be < 1.4 for random configurations. 

It is the purpose of this note to show that (C) can be 
easily estimated by using well known formulae for generic 
intersections of random fractals |Tc|] . Take two fractal 
sets X and Y with dimensions Dx and Dy, embedded 
in a space of dimension d. Then, their intersection has 
dimension DxnY = Dx + Dy — d for nearly every rel- 
ative position and orientation, provided this intersection 
is non-empty. In the present case, X is the curve to be 
studied, Y is a line of view (therefore Dy = 1) which 
passes through X, and d = 3. This gives 



D 



xnY 



= D x -2. 



(2) 



If this is positive, the average number of intersections 
between the line of view and X increases as m ~ 
N D X nY/D x for — s- oo. Actually, for this to be true 
we either have to assume that X is a true fractal without 
lower length cutoff (which is not true for random walks 
with finite step size a), or we have to fatten Y. Thus we 
consider instead of a single line of view a cylinder whose 
thickness is of the order of the step size a (the precise 
value is irrelevant), and the above number of intersec- 
tions has to be interpreted as the number of crossings 
between projected bonds within a distance 0(a). This 
scales in the same way as the number of crossings per 
bond. Thus we obtain immediately 



with 



{C)/N~N & (3) 



1 = Dxnv/Dx = 1 - 2/D x . (4) 



This is the case for compact (collapsed) polymers where 
Dx = 3 and therefore 



a = 1/3 (compact globules), 



(5) 



in perfect agreement with the exact results of |ll| , P-2| and 
with the numerical values for protein backbones. 
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This argument has to be modified when Dx < 2. 
For D x < 2 (i.e., for SAWs where D S aw = 1/v with 
v = 0.5877 ]13|]), the above argument gives, mutatis 
mutandis, not the leading term but the first sublead- 
ing correction. The leading term is (C)/N ~ const 
since {C)/N cannot vanish for N — > oo The av- 

erage number of intersections per bond scales now as 
to ~ const - N DxnY / Dx . This gives 

(<7)/iV ~ a - ai/AT (SAW) (6) 

with a given by Eq.(^). 

For Dx = 2, finally, we expect logarithmic depen- 
dence. This is the case for ordinary random walks (RW) 
and for 9-polymers (actually, since we want no true self- 
intersections, we should consider here only the latter. 
But we shall speak of RWs for simplicity. Later, in eq.(||), 
we will also allow true RWs). The number of times a 
projection of an TV-step RW comes back to a previously 
visited site, within a finite distance a, increases ~ In AT 
Q , and therefore 

(C)/N~hiN. (RW) (7) 

In order to verify Eqs.^J^) numerically, we performed 
Monte Carlo simulations (we did not make simulations 
for the collapsed case since there the theoretical result is 
too obvious). In these simulations we indeed did not cal- 
culate (C) proper, but a closely related quantity which 
should show the same asymptotic behaviour and is much 
easier to calculate numerically. 

We consider walks on a simple cubic lattice, and con- 
sider only projections along one of the 3 axes. We con- 
sider the sites i on the projection plane which arc visited 
irii times with m, > 1, and count the number of pairs of 
visits to the same site. Dividing this by N + 1 , we obtain 

^^E^-^ (8) 

where i runs over all sites the plane. We call this the 
"opacity" . It can obviously also be defined for ordinary 
(i.e. non-self avoiding) walks. If there are no double vis- 
its in the projection, the opacity is zero, while it diverges 
with N for a compact object whose thickness along the 
line of view diverges with N . 

The main differences of B with respect to (C) are 
that we do not require transversality of the crossings and 
do not make any angular averaging. Indeed we project 
along atypical directions where bonds would not intersect 
transversally but overlap, so that C cannot be properly 
defined. As a consequence, the numerical calculation of 
B is trivial in comparison with the calculation of (C), 
for any given configuration. Nevertheless, we conjecture 
that both show the same asymptotic behaviour. Indeed, 
if this were not the case (i.e., if C would depend strongly 



on the angle of view), (C) would be an average over a 
strongly fluctuating quantity and would presumably not 
be of much practical use. But there exists strong evi- 
dence that projections and intersections of random frac- 
tal objects show generic features independent of the angle 
of projection resp. intersection. Finally, our theoretical 
discussion used generic projections and applies therefore 
strictly spoken only to (C). If we find agreement for B, 
this suggests that the arguments are correct for (C) a 
fortiori. 




N 



FIG. 1. Log-linear plot of B versus N for ordinary ran- 
dom walks (full line). The dashed line shows the analytic 
prediction Eq.(^), with an additive constant not predicted 
by theory. The value of this constant obtained by fitting the 
data is -0.398 ± 0.03. 
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FIG. 2. Log-linear plot of B versus N for SAWs (full 
line). The dashed line shows the analytic prediction Eq.(^), 
where we have fitted numerically the constants ao — 1.39 
and a± — 1.415. 
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For RW, we can indeed estimate B more precisely, 
and give the prefactor in the analogon to eq. ([?]). A pro- 
jection of an iV-step lattice RW along one of the coordi- 
nate axes is a RW of 2JV/3 steps in the plane, for large 
N . The number of distinct sites visited by the latter is 
w (2wN/Z)/hiN @. Thus m « (3/27r)ln7V, and 

B~-£-hiN. (RW) (9) 

Results from 10 4 RWs of 4 x 10 5 steps each are shown 
in Fig.l, together with the prediction Eq.(||). Adding an 
offset which we had not tried to calculate analytically 
and which turned out to be precisely —0.4 within the 
estimated error bars, we find perfect agreement. 

For simulating SAWs we used the PERM algorithm 
p5[ with Markovian anticipation bias |lq ]. This sam- 
ple contained 5 x 10 6 SAWs (of which ca. 150,000 
were strictly independent) of length 4000. Results are 
shown in Fig. 2. There we show also the curve 1.39 — 
1.415/iV 54 which obviously gives a perfect fit for large 
N, verifying Eq.(^). We should point out that our numer- 
ical values are in surprisingly close agreement with those 
shown in Fig. 3 of given the fact that we do not mea- 
sure exactly the same quantity. In particular, if we would 
make a least square fit to our data with 400 < N < 1500, 
we would also get a ~ 1.1, in rough agreement with [Q. 
But clearly such a fit would have a disastrous chi-squared. 
This clearly suggests that corrections to normal scaling 
have been mistaken in Jl]|| for an anomalous power law. 

Up to now we have only considered walks on the sim- 
ple cubic lattice. Off-lattice polymers can be treated in 
the same way, by replacing the condition of exact coin- 
cidence by an approximate one. Two monomers i and 
k contribute then to B(n), where h is the direction of 
projection, if \{xi — X]~)h\ < e for some suitably chosen 
accuracy e. Notice that the scaling behaviour of B{n) 
should not depend on e. We should add finally that B 
can also be evaluated for any set of sites, not necessar- 
ily being lined up to form a toplogically linear chain. 
We can use it therefore to measure opacities of branched 
polymers, vesicles, clusters, or droplets. 

Although we have not proven rigorously that B scales 
in the same way as (C) /N, we believe that there is little 
room to doubt it. In any case, neither (C) nor B are 
proper measures of entanglement, since e.g. the same 
scaling laws as for open chains should be observed also 
for unknotted loops. Instead, both are measures of opac- 
ity. Of course they can be used to monitor the theta 



transition from an open coil to a collapsed globule, since 
the transparency of a coil decreases during the collapse. 
But it is not clear what advantage they offer compared, 
e.g., to the gyration radius. The large corrections to the 
asymptotic behaviour seen in Fig. 2 (which have previ- 
ous authors even mislead to postulate anomalous scal- 
ing laws) should be a warning that the interpretation of 
numerical values might be difficult sometimes. In any 
case, since the calculation of B is much simpler than 
that of (C) while it gives basically the same information, 
it should be prefered in any application. 

I am indebted to Walter Nadler and Rodrigo Quian 
Quiroga for very useful and stimulating discussions, and 
to Stuart Whittington for pointing out an error in the 
original manuscript. 
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